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Abstract 

Energy harvesting (EH) provides a means of greatly enhancing the lifetime of wireless sensor 
nodes. However, the randomness inherent in the EH process may cause significant delay for performing 
sensing operation and transmitting the sensed information to the sink. Unlike most existing studies on 
the delay performance of EH sensor networks, where only the energy consumption of transmission is 
considered, we consider the energy costs of both sensing and transmission. Specifically, we consider an 
EH sensor that monitors some status property and adopts a harvest-then-use protocol to perform sensing 
and transmission. To comprehensively study the delay performance, we consider two complementary 
metrics and analytically derive their statistics: (i) update age - measuring the time taken from when 
information is obtained by the sensor to when the sensed information is successfully transmitted to the 
sink, i.e., how timely the updated information at the sink is, and (ii) update cycle - measuring the time 
duration between two consecutive successful transmissions, i.e., how frequently the information at the 
sink is updated. Our results show that the consideration of sensing energy cost leads to an important 
tradeoff between the two metrics: more frequent updates result in less timely information available at 
the sink. 
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transmission. 
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I. Introduction 

Background: Energy harvesting (EH) from energy sourees in the ambient environment is an 
attraetive solution to power wireless sensor networks (WSNs). The feasibility of powering WSNs 
by EH from solar, wind, vibration and radio-frequeney (RE) signals has been demonstrated in 
the literature [[T]-[^. If an EH souree is periodieally or eontinuously available, a sensor node ean 
in theory be powered perpetually. However, the design of EH WSNs raises several interesting 
and ehallenging issues. 

Design Challenges: An important design eonsideration for EH WSNs is the modeling of 
energy eosts. There are three main energy eosts in wireless sensors Q: (i) energy eost of 
RE transmission and reception, including idle listening, (ii) energy cost of information sensing 
and processing, and (iii) energy cost of other basic processing while being active. Generally, 
the energy cost of other basic processing is much smaller compared to the energy cost of 
transmission 0 - Hence, the majority of the current work on EH WSNs has considered only 
the energy cost of transmission, while ignoring the energy cost of sensing Q, Q. Eor some 
sensors, such as high-rate and high-resolution acoustic and seismic sensors, the energy cost of 


sensing can actually be higher than the energy cost of transmission, e.g., see [ 101 and references 
there in. Hence, it is important to accurately model the energy cost of sensing in WSNs pT| . 

Eor WSNs powered by EH from the ambient environment, the energy arrival process is 
inherently time-varying in nature. These fluctuations in the energy arrival process can be slow 


or fast and are characterised by its coherence time [12|. Eor instance, for the case of EH from a 
solar panel on a clear day with abundant sunshine, the coherence time is on the order of minutes 
or hours. Eor the case of wireless energy transfer via RE signals, the coherence time can be on 
the order of milliseconds, which is comparable to the duration of a communication time slot. 
The energy arrival process in the latter case can be modeled as a random process where the 
amount of harvested energy in each time slot follows some probability distribution. Eor example, 
papers studying EH from RE signals often assume an exponential distribution [13 -1^. Another 
example, using the gamma distribution, can be found in [ fT^ . However, many energy arrival 
processes in practice cannot be accurately modeled by using exponential or gamma distributions. 
The consideration of a more general probability distribution for modeling the amount of energy 
arrival is still largely an open problem. 
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In many sensor network applieations, the delay performanee is a key design ehallenge. The 
etfeets of randomness in both arrivals of the multiple data paekets and harvested energy on 


the overall transmission eompletion time were eonsidered in [17|. In a single data paeket 
and randomness in the energy arrival proeess and wireless ehannel, were eonsidered in the 
analysis of transmission delay, i.e., the time duration between the generation of a paeket and its 
sueeessful reeeption at the sink. However, both and [ 13| only eonsidered the energy eost of 
transmission. To the best of our knowledge, a eomprehensive analysis of the delay performanee of 
EH WSNs taking into aeeount a realistie model of sensor energy eosts, has not been investigated 
in the literature. 


Paper Contributions: We eonsider a status monitoring seenario, e.g., monitoring some prop¬ 
erty of a target environment, with one sensor-sink pair. The sensor is solely powered by EH from 
an ambient energy souree. The sensor periodieally monitors and senses the eurrent environment, 
i.e., it generates eurrent status information about one or more variables of interest, and then 
transmits a status-information-eontaining paeket to the sinl|^ Onee the paekets are sueeessfully 
transmitted to the sink, whieh may oeeur after several failed retransmissions due to fading in 
the transmission ehannel, the status under monitoring is updated at the sink. 

We adopt two different metries to assess the delay performanee: (i) update agej^whieh measures 
the time duration between the time of generation of the eurrent status information at the sink 
and the time at whieh it is updated at the sink, and (ii) update eyele whieh measures the time 
duration between one status update at the sink to the next. The update age (or freshness) and 
update eyele (or frequeney) are eomplementary measures. Eor instanee, a smaller update age 
means the updated status information at the sink is mueh more timely, but does not indieate when 
the next update status information will be reeeived. A smaller update eyele means more frequent 
status updates at the sink, but does not indieate when the eurrent updated status information 


'Due to the fluctuation in the energy arrival process, strictly periodic sensing and transmission is not possible. In this paper, 
‘periodic’ is used to indicate that the sensor alternates between sensing and transmission(s) in order to keep status updating at 
the sink. 


The term update age is inspired by |18| and indicates the age or timeliness of the transmitted information, since an outdated 
message may lose its value in a communication system when the receiver has interest in fresh information (T9) . Note that this 
notion of the delay is in fact the same as the transmission delay in |13|. 
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was originally generated or how old it is. Thus, the quality of a status monitoring system, i.e., 
the status update freshness and frequeney, is eomprehensively eaptured by the update age and 
update eyele, respeetively. 

We aeeount for the faet that sensing and transmission operations both eonsume energy. Inspired 
from the harvest-then-use and save-then-transmit eommunieation protoeols for EH nodes in 
wireless networks p3| , [14|, [16|, whieh are simple to implement in praetiee, we eonsider 
a harvest-then-use protoeol for the EH sensor. In our proposed protoeol, the sensor performs 
sensing and transmission as soon as it has harvested suffieient energy. In order to limit the delay 
due to retransmissions, we impose a time window for retransmissions. The delay performanee 
of the eonsidered harvest-then-use protoeol is analyzed. The main eontributions of this paper are 
as follows: 

• We provide a eomprehensive study on the delay performanee of EH sensor networks. Apart 
from the eommonly eonsidered delay due to the information transmission from the sensor 
to the sink, defined as the update age, we also eharaeterize the frequeney of updating the 
information held by the sink, defined as the update eyele. 

• Considering a Rayleigh fading wireless ehannel, we analytieally derive the statisties of both 
the update eyele and the update age. We eonsider both a deterministie energy arrival model 
and a random energy arrival model with a general distribution, so that our results ean be 
applied to model a wide range of EH proeesses. 

• We take the energy eosts of both sensing and transmission into aeeount when studying the 
delay performanee. Sueh a eonsideration brings up an interesting question of whether to in- 
erease or reduee the number of allowed retransmission attempts for eaeh sensed information, 
beeause both sensing and transmission eonsume energy. This in turn results in a tradeoff 
between the update eyele and the update age. The tradeoff emphasizes the importanee of 
modeling the energy eost of sensing. 

Notations: IE{-} and Pr{-} are expeetation and probability operators, respeetively. Convolu¬ 
tion operators for eontinuous and diserete funetions are denoted as and *, respeetively. [■] and 
[•J are eeiling and floor operators, respeetively. summation operator, and if m > n, 

the result is zero. Pois {i, A) is the probability mass funetion (pmf) of a Poisson distribution with 
parameter A. 




5 


II. System Model 

We consider the transmission scenario where a sensor periodically transmits its sensed in¬ 
formation to a sink, as illustrated in Fig. The sensor is an EH node which harvests energy 
from the ambient environment such as solar, wind, vibration or RF signals. The sensor has two 
main functions, i.e., sensing and transmission, each having individual energy cost. We assume 
half-duplex operation, i.e., sensing and transmission cannot occur at the same time. In order to 
perform either sensing or transmission, the sensor first needs to spend a certain amount of time 
on EH. The harvested energy is stored in a battery. We assume that the battery cannot charge and 


discharge at the same time [16|. In addition, the battery has sufficient charge capacity such that 
the amount of energy stored in the battery never reaches its maximum capacity. This assumption 
is reasonable since battery capacity typically ranges from joules to thousands of joules [|T||, while 
the energy level in the battery in our system is only in the /iJ range as shown in Section V. 


Following the state-of-the-art EH sensor design practice [20|, we adopt a time-slotted or block- 
wise operation. We assume that one sensing operation or one transmission is performed in one 
time block of duration T seconds]^ At the beginning of each block, we assume that the sensor 
checks the battery energy state and makes a decision to perform either sensing, transmission, 
or energy harvesting. Thus, we define the following types of time blocks with the associated 
amount of energy cost/harvesting: 

• Sensing Block (SB): the sensor samples the status information and then processes and packs 
sensed information into a data packet. The energy cost in a SB is denoted by £^sb- 

• Transmission Block (TB): the sensor transmits the newest generated data packet (from the last 
sensing operation) to the sink with energy cost i-C-, the transmit power is Ptb = ^tb/T. 
Then the sink sends a one-bit feedback signal to the sensor to indicate successful packet 
reception. We assume that the time consumed for receiving the feedback signal at the sensor 
is negligible as compared to its packet transmission time. If the transmission is successful, we 
have a successful transmission block (STB); otherwise, we have a failed transmission block 


(FTB). We assume that successes/failures of each TB are mutually independent [13|, [14|. 
The probability of a TB being a FTB, i.e., transmission outage, is denoted by Pont- 


^In general a sensor may spend different amounts of time on one sensing operation jlo[ . Thus, the assumed protocol and 
analysis can be generalized to different sensing time durations other than T, which is outside the scope of this work. 
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Fig. 1 Illustration of system model and sensor components. 

Energy-harvesting bloek (EHB): the sensor harvests energy from the ambient environment and 
stores the energy in its battery. 


A. Proposed Sensing and Transmission Protocol 

Since the time-varying EH process results in randomness in the delay for performing sensing 
and transmission, we propose a harvest-then-use protocol with a time window for retransmissions 
in order to improve the delay-related performance. 

The protocol is motivated as follows. Firstly, considering the energy cost of sensing, it is 
necessary to harvest sufficient energy, before sensing can occur. However, it is unwise to 
perform sensing as soon as the harvested energy reaches £^sb because there will be insufficient 
energy left for transmission after the sensing operation. The time spent on EH due to insufficient 
energy for transmitting the sensed information will result in unnecessary delay. To avoid such 
delay, we define the condition for the sensing operation to be when the harvested energy in the 
battery exceeds £^sB + i^TB- In this way, a transmission of sensed information occurs immediately 
after the sensing operation (i.e., a SB is always followed by a TB). Secondly, in the event 
that the transmission is not successful due to the fading channel between the sensor and sink, 
we need to allow for retransmissions, which are a common feature in conventional (non-EH) 


WSNs [211. In this paper, we impose a time window for retransmissions to control the delay 


caused by unsuccessful transmissions because it is unwise to spend an indefinite amount of time 
trying to transmit outdated information. We denote W as the maximum number of time blocks 
after a SB, within which transmissions of the currently sensed information can take place. Since 
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Fig. 2 Illustration of update cycle and update age. 


the first transmission attempt always happens immediately after the SB, the time window for 
retransmissions is fU — 1 time blocks. 

Under the proposed protocol, the sensor operates as follows: 

1) First, the sensor uses several FHBs to harvest enough energy, £^sb + and then a SB 
and a TB occur. 

2) If the transmission in the TB is successful, i.e., we have a STB, the sensor harvests energy 
(taking several FHBs) for the next sensing period until the battery energy exceeds £^sb+^tb- 

3) If the transmission in the TB fails, i.e., we have a FTB, the sensor goes back to harvest 
energy (taking several FHBs) and performs a retransmission when the battery energy 
exceeds £^tb- 

4) Retransmission may occur several times until the sensed information is successfully trans¬ 
mitted to the sink or the time window for retransmissions W — 1 is reached. Then, the data 
packet at the sensor is dropped and the sensor goes back to harvest the energy for a new 
sensing operation. 

Fig. illustrates this protocol with IF = 7. In the example shown, the first block in Fig. 
is a SB, followed by two FTBs (and two FHBs in between). Since the third TB is a STB, the 
sensed information in the first SB is successfully transmitted to the sink. Then, the sensor uses 
three FHBs to harvest energy to conduct sensing in the next SB. After the second SB, there 
are three TBs during 7 time blocks, and all of them are FTBs. Thus, the retransmission process 
is terminated after IF = 7 is reached. As a result, the sensed information in the second SB is 
not transmitted to the sink. The time indices shown in Fig. will be defined in the following 
















section. 


B. Proposed Models for Energy Arrival 

In this paper, we consider that the harvested energy in each EHB could either remain constant 
or change from block to block. The former is referred to as deterministic energy arrival, while 
the latter is referred to as random energy arrival. 

Deterministic energy arrival is an appropriate model when the coherence time of the EH 
process is much larger than the duration of the entire communication session, such as EH by 


solar panel on clear days p^ , [ |22| , p3| |. In this paper, we denote this as deterministic energy 
arrival process. Eor tractability, we also assume that £’sb and £tb represent integer multiples of 
the harvested energy by one EHB, p. 

Eor random energy arrivals, we consider independent and identically distributed (i.i.d.) random 
energy arrival modeQ with a general probability distribution function for the amount of energy 
harvested in each EHB. This energy arrival model is referred to as general random energy 


arrival process. The previously considered exponential and gamma distributions in [13|, [14|, 
CD become as special cases of the general probability distribution in this work. Since the 
exponential distribution is commonly studied for wireless power transfer using RE signals, we 
will also provide results for this important special case and referred to it as exponential energy 
arrival process. 


HI. Delay-Related Metrics 

As described in the previous section, both sensing and (re)transmission requires a variable 
amount of EH time, which may result in significant delays in obtaining the sensed information 
at the sink. In this section, we consider two metrics to measure the delay performance of the 
considered sensing and transmission protocol. 

Eor the convenience of describing the two metrics, as shown in Eig. we use fsTBj to denote 
the block index for the jth STB during the entire sensing and transmission operation. Note that 
a successful transmission also induces an information update at the sink. Also, it is important 


''The i.i.d. energy arrival model is commonly considered in the literature 


0 . 


25 


. There are other energy arrival models 


captures the temporal correlation of the energy arrival process, such as discrete-Markovian modeling Q, (0, |26| , which are 
beyond the scope of this work. 
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to associate each transmission with its information content. To this end, we use tsBj to denote 
the block index for the SB in which the sensed information is transmitted in the jth STB. In 
other words, status information sensed at fsBj is successfully transmitted to the sink at fsTBj- 
Next, we define two delay-related metrics, expressed in terms of the number of time blocks: 

A. Update Age and Update Cycle 

Definition 1 (Update age). For the jth STB, the update age is given by the number of time 
blocks from fsBy to fsTBj (shown in Fig. [^. The jth update age is 

^UAj = tsTBJ — tsB,j, j = 1, 2, 3, .... (1) 

Remark 1. The update age measures the time elapsed from the generation of a status-information- 
containing packet at the sensor to the reception of the packet, i.e., status update, at the sink. 
This metric is referred to as the status update age in p8| . A larger update age implies that a 
more outdated status is received by the sink. The update age, which captures the freshness of the 
updated status information, however, does not reflect the update frequency at the sink. Rather, 
the update frequency is captured by the update cycle which is presented below: 

Definition 2 (Update cycle). For the jth STB, the update cycle is given by the number of time 
blocks from fsTBj-i to fsTBj (shown in Fig. [^. The jth update cycle is 

Tvc,j = tsTB,j+l — tsTB,j, j = 1, 2, 3, .... (2) 

Remark 2. The update cycle measures the time elapsed from one status update at the sink to 
the next. The update cycle, however, does not reflect the update freshness at the sink. Unlike 
the update age, the update cycle takes into account the delay due to dropped data packets. 
Therefore, update cycle complements update age, and they jointly capture the update frequency 
and freshness, to provide comprehensive metrics on the delay performance of a status monitoring 
system. 

B. Modeling Delay-Related Metrics as i.i.d. Random Variables 

To model each of the update age/update cycle as i.i.d. random variables, we focus on the 
steady-state behavior as characterized in Lemma 1. 
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Lemma 1. For a deterministic energy arrival process, the energy level after each TB is zero. For 
a general random energy arrival process with pdf containing at least one positive right-continuous 
point, /(e), the steady-state distribution of the energy level after each TB has pdf 

9(£) = i(l-F(£)), (3) 

where p is the average harvested energy, and F{e) is the cumulative distribution function (cdf) 
corresponding to /(e). 

Proof: For a deterministic energy arrival process. Lemma 1 is straightforward. For a general 
random energy arrival process, the proof is given in Appendix B. ■ 

According to the sensing and transmission protocol defined in the previous section, each SB 
is directly followed by a TB. From Lemma 1, the steady-state distribution of available energy 
after any TB is the same. Hence, the steady-state distribution of the available energy after fsTBj 
is the same for all j. Because the successes/failures of each TB are mutually independent, and 
Tucj is determined by both the available energy after fsTBj and the successes/failures of the 
following TBs, Tucy are i.i.d. for all j. Similarly, it is also easy to show that Tuaj are i.i.d. 
for all j. For convenience, we remove subscript j for Tuc and Tua in @ and ([T]), respectively. 

IV. Update Age 

In this section, considering the dynamics of an energy arrival process and the probability 
of successful/failed transmission in our proposed harvest-then-use protocol, the update age for 
deterministic, general random and exponential energy arrival processes are analyzed. 


A. Deterministic Energy Arrival Process 

Theorem 1. For a deterministic energy arrival process, the update age pmf is given by 

£tb 


ID 1 ^ -Pout) (-Pout) J 1 I /' 

Pr{TuA = K|=---, k = l + [n-l) 


where 


Ps, 


n = 1, 2, ...h, fi = 1 + 


P 


+ 1 , 


lU- 1 

^TB 

. ' P . 


, -Psuc = 1 — (-Pc 


out ) ? 


and Pout is the probability of a TB being a FTB, defined in Section II. 


(4) 

(5) 


Proof: See Appendix D. 
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From Theorem 1, the average update age, Tua for a deterministie energy arrival proeess is 
straightforwardly obtained as in Corollary 1. 

Corollary 1. For a deterministic energy arrival process, average update age is given by 

ruA - E (i+(" -1) (^+1)) ~ 

where Pguc is given in 


B. General Random Energy Arrival Process 

Theorem 2. For a general random energy arrival process, the update age pm^ is given by 
ri-Pout 


P., 


k = 1, 


Pr{TuA=fc}=< 


where 


(1 Pout)^ - l)^TB)-G'fe-n((n - 1)^tb)) ,2<k<W, 

(7) 


P., 


n=2 


I 


P,,c = 1 - Pout + (1 - Pout) (Pout)"-' - 1)^tb) - Gi.ni{n - 1)^tb)) , 


1=2 n=2 


and 


GAx) = 


i = -1, 


( 8 ) 


(9) 


V. 0 i convolutions 

g{x) and f{x) are defined in Lemma 1. 

Proof: See Appendix D. ■ 

From Theorem 2, the average update age, Tua for a general random energy arrival proeess is 
obtained straightforwardly as in Corollary 2. 

Corollary 2. For a general random energy arrival process, average update age is given by 

w i 


Tua = 


1-P 


out 


P, 


l + 5^/5^(Pout)”"'(G';-n-l((n-l)^TB)-G'l-n((n-l)^TB)) , (10) 


1=2 n=2 


^Although the general expression in Theorem 2 contains multiple integrals in Eq. (9), for special cases, such as deterministic 
and exponential energy arrival processeses, the results given in Theorems 1 and 3 are closed-form expressions. 
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where Psuc is given in 


C. Exponential Energy Arrival Process 


Theorem 3. Eor an exponential energy arrival process, the update age pnif is given by 
( 1 - Pont 


Pr {Tua = k} = < 


(1-^c 


-,/c = 1, 


out J 


^ Pois {k -n,{n- 1)£:tb/p) ,“2 <k<W, 


( 11 ) 


n=2 


where 


w i 


Psuc — 1 ~ -Pout + (1 ~ -Pout) EE (-Pout)" ^ Pois(/- n, (n - 1 )£:tb/p) , (12) 


1=2 n=2 


Proof: See Appendix D. ■ 

From Theorem 3, the average update age, Tua for an exponential energy arrival process is 
straightforwardly obtained as in Corollary 3. 


Corollary 3. Eor an exponential energy arrival process, average update age is given by 

tuA = \ ^ + {PoutT~" Pois {I -n,{n- 1)Ptb/p) ] (13) 

V 1=2 n =2 J 

where Psuc is given in & 

From Theorems 1-3 and Corollaries 1-3, we see that different energy arrival processes induce 
different pmfs and average values of update age. For benchmarking with the existing studies 
on delay without imposing a constraint on the time window for retransmissions pA| , we let 
W —)■ cx), so that all sensed information is eventually transmitted to the sink, the average update 
age is the same under different energy arrival processes as in Corollary 4. 

Corollary 4. Eor a deterministic or general random energy arrival process, Tja increases with 
W, and as W gets large, the asymptotic upper bound of Tua is independent with energy arrival 
distribution and is given by 

hm Tua = 1 + f— + iV (14) 

W^oo 1 - Pom \ p J 


Proof: See Appendix G. 
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Remark 3. From the above analytieal results, we have that: 

i) From Theorems 1 to 3, Tua is independent of the energy eost of sensing, beeause 
the delay is only affeeted by the energy harvesting and retransmissions that happen after 
the sensing operation. This might give the impression that energy eost of sensing does not 
affeet delay. However, update age is only one of the two delay metries, and the energy eost 
of sensing has important impaets on update eyele, whieh will be investigated in the next 
seetion. 

ii) Allowing a larger window for retransmissions inereases the average update age. This might 
suggest that retransmissions should be avoided, i.e., W = 1. However, the update age does 
not take into aeeount oases where sensed information is not suooessfully transmitted to the 
sink. In this regard, the update eyele implioitly oaptures sueh oases. 

V. Update Cycle 

In this seetion, oonsidering the dynamios of an energy arrival proeess and the probability of 
sueeessful/failed transmission in our proposed harvest-then-use protoeol, the update eyele for 
deterministio, general random and exponential energy arrival prooesses are analyzed. 


A. Deterministic Energy Arrival Process 

Theorem 4. For a deterministic energy arrival process, the update cycle pmf is given by 

Pr {Ac = ^ = (1 - A.0 (Pou.)-"’”*. k = + (ft + 1 )| 

(15) 

where n = 1, 2, ...n, m = 0, 1,2,..., and n is given in Q. 

Proof: See Appendix E. ■ 


Corollary 5. For a deterministic energy arrival process, average of update cycle is given by 
(Pout)" 


Puc = 


(Pc 


out ) 


1 + ft + + Ai +1+(1 + JaA 

P J P P l-(Pc 


out 

h 

out; n=l 


(Pout)" ^n. 

(16) 


Proof: See Appendix F. 
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B. General Random Energy Arrival Process 

Theorem 5. For a general random energy arrival process, the update cycle pmf is given by 
Pt {Tuc = k} =y^ic{EsB+£TB) *C{£sb) * • • • * C(^sb)^*^ * V ' * + k = 2,3, 

m=o\ m convolutions m convolutions J 

(17) 

where rh = and functions C,{£,i), i{i) and -d^i) are given by 

a£,t) = G,.i{£)-Gi{£), (18a) 

L{i) = PsucPr {TvA = i} , (18b) 

'd{i) = Pout {Gw+i-2{£-TB) — Gw+i-i{£tb)) 

W I 

+ 'y~^'y~^ (Pout)” {Gl-n-l{{n — 1 )£tb) —Gi-n{{n — 1)£^tb)) {Gw+i-l-l{^TB) —Gw+i-l{£TB)) ■ 

1=2 n=2 

(18c) 

Pr{TuA = i}, Psuc cind Gi{£) are given in Q, ® and respectively. 

Proof: See Appendix E. ■ 


Corollary 6. For a general random energy arrival process, average update cycle is given by 

fvc = ^ ~ (— + V + W + l] + ^ + Tua + 1, (19) 

Psuc VP J P 

where Psuc and Tua are respectively given in @ and eg, and 

f W-2 / W-2 

V = T——— I- i (Gj_i(£^TB) — ^^(f^TB)) — (IP— 1) 11— (Gj-i (£^tb) — Gi{£-rB)) 

\ P i=0 \ i=0 


w i 


1-P, 


(Pout)"^ {Gl-n-l{{n — l)i^TB) — Gi_n{{n — 1)^^tb)) 


X 


1=2 n=2 


& 


W-l-l 


W-l-1 


TB 


P 


* (G',_i(^tb) - ^.(^tb)) - (IP - /) 1- (G',_i(^tb) - Gi{£TB)) 


i=0 


i=0 


( 20 ) 


Proof: See Appendix F. 
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C. Exponential Energy Arrival Process 

Theorem 6. Eor an exponential energy arrival process, the update cycle pmf is given by 

— m(l + W) — 1), k = 2,3,.... 

(21) 

where m = and functions C,{S,i), i{i) and -§{ 1 ) are given by 

C,{S,i) =Poi?,{i,8/p), (22a) 

i{i) = -PsucPr {Tua = i} , (22b) 

d{i) = PoutPois {W + i- l, Stb/p) 

w i (22c) 

(Pout)” Pols (/ -n,{n- 1)£^tb/p) Pols (W + i - I, Stb/p) , 

1=2 n=2 

and Pr{TuA = i} cind Pguc are given in and ( [T^ , respectively. 

Proof: See Appendix E. ■ 


m / 

Pr {Tuc = k} = ^ ( C(("i + 1)^^SB + ^^tb) if .y if 

771=0 \ m convolutions 


Corollary 7. Eor an exponential energy arrival process, average update cycle is given by 


f(^ + + iv + 1 ) + + f „* +1, 

Psuc VP J P 


(23) 


where Tua cind Pguc cire given in and iB. and 

W-2 


P = 


P 


out 




TB 


1-P., 


P 


i=0 


zPois ( i, j — {W — 1) I 1 — ^ Pois ( i. 


P 


w -2 


i=0 


P 


TB 

P 


1-P, 


-X 


W I 


(Pout)” Pois ( l-n, (n-1) 


1=2 n=2 


P 


TB 


P 


'P 


W-l-1 


TB 


P 


i=0 




P 


W-l-l 


^^iPois ( i, j —(IP —Z) |l— ^Pois ( i. 


i=0 


P 


TB 


(24) 


Proof: See Appendix F. ■ 

Similar with the case of update age, different energy arrival proeesses induees different pmfs 
and average values of update eyele. However, for benchmarking with the existing studies on 
delay without imposing a eonstraint on the maximum allowable retransmission time, when we 
consider removing the constraint of retransmission, i.e., IP —)• 00 , so that all sensed information 
is eventually transmitted to the sink, the average update cycle is the same under different energy 
arrival proeesses as in Corollary 8. 
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Corollary 8. For a deterministic or general random energy arrival process, T\jq decreases with 
W, and as W grows large, the asymptotic lower bound o/Tuc is independent with energy arrival 
distribution and is given by 

y rf, r, , ^^SB + ^TB , -Pout /Ptb , 

hm Tuc = 2 +-+ --— -+ 1 . (25) 

W^oo p 1 - \ p J 

Proof: See Appendix G. ■ 


Remark 4. From the above analytieal results, we have that: 

i) From Theorems 4 to 6, we know that Tuc is affeeted by the energy eost of sensing, £^sb- 
A larger £^sb means more EHBs are required to harvest a suffieient amount of energy to 
perform sensing operation(s) between adjaeent STBs. 

ii) A larger window for retransmissions shorten the average update eyele, beeause allowing 
more retransmissions inereases the ehanee of having a sueeessful transmission. This might 
suggest that it is also better to inerease W to reduee the update eyele. But inereasing W also 
inereases the update age as diseussed earlier. Therefore, there is elearly a tradeoff between 
the two metries. 


VI. Numerical Results 

In this seetion, we present numerieal results for the update age and update eyele, using the 
results in Theorems 1-6 and Corollaries 1-8. The typieal outdoor range for a wireless sensor is 


from 75 m to 100 m [27|. Henee, we set the distanee between the sensor and the sink as d = 90 m 


and the path loss exponent for the sensor-sink transmission link as A = 3 [ 13|. The duration of a 
time bloek is T = 5 ms [22|. The noise power at the sink is cx^ = —100 dBm [14|. The average 


harvested power is 10 mW [ [28| , i.e., average harvested energy per time bloek, p = 50 pi. 
Unless otherwise stated, (i) we set the power eonsumption in eaeh TB, Ptb = 40 mW, i.e., 
Ptb = 200 pi. Note that this ineludes RF eireuit eonsumption (main eonsumption) and the 
aetual RF transmit power Vt^ = —5 dBmj^and (ii) we set the power eonsumption in eaeh SB as 


PsB = 50 mW p0| , i.e., Psb = 250 pi. In the following ealeulations, power and SNR related 
quantities use a linear seale. We assume that a transmission outage from the sensor to the sink 


"The values we chose for Ptb and Ptx are typical for commercial sensor platforms, such as MICAz 


0 
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occurs when the SNR at the sink 7 , is lower than SNR threshold 70 
probability is 

Pout = Pr {7 < 70 } . 


40 dB p5| . The outage 
(26) 


The SNR at the sink is |29| 


7 = 




tx 


rdV2 


(27) 


where h is the souree-sink ehannel fading gain, P = is a path loss faetor relative to 

referenee distanee do of the antenna far field, and Piido) is linear-seale path loss, whieh depends 
on the propagation environment Following p3| , pl| , we assume P = 1, for simplieity. 

For the numerieal results, we assume that h is bloek-wise Rayleigh fading. Using ( [27] ), the 
outage probability ean be written as 


P 


out 


1 — exp 




(28) 


By applying ( [28| ) to the theorems and eorollaries in Seetions IV and V, we eompute the expres¬ 
sions of the pmfs of Tua and Tuc as well as their average values 7 \ja and 7\jc- 

Pmfs of update age and update cycle with different energy arrival processes: First, we 
eonsider a deterministie energy arrival proeess with harvested energy in eaeh EHB, p. Also 
we eonsider two speeial oases of the general random energy arrival proeess: (i) exponential 
energy arrival proeesses with average harvested energy in eaeh EHB, p and (ii) random energy 
arrival proeesses with gamma distribution [p^ , Gamma(0.05,1000). We term this as the gamma 
energy arrival process, and it is easy to verify that this gamma energy arrival proeess has the 
same average harvested energy in eaeh EHB as the deterministie and exponential energy arrival 
proeesses. 

Figs. |3]j^ plot the pmfs of update age, Tua, and update oyole, Tuc, for the deterministie, 
gamma and exponential energy arrival proeess, respeetively. The analytieal results are plotted 
using Theorems 1-6, and we set W = 50, i.e., the time window for retransmissions is lU — 1 = 
49 time blooks. In partioular, in Fig. the analytieal pmfs of Tua and Tuc for the general 
random arrival proeess are obtained using Theorems 2 and 5. The results in Figs. |4]j^ also 
illustrate the importanee of the general random energy arrival proeess, whieh is used in this 
work. This is beeause gamma and exponential energy arrival proeesses, whieh have been used in 


the literature [ 13 -161, are speeial eases of the general random energy arrival proeess. We see that 
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3 60 80 100 120 140 160 180 

T\jq, time blocks 


Fig. 3: pmfs for Tua and Tuc with Fig. 4: pmfs for Tua and Tuc with Fig. 5: pmfs for Tua and Tyc with 
deterministic energy arrival process, gamma energy arrival process. exponential energy aiTival process. 


different energy arrival processes result in different pmfs of update age and update cycle. Hence, 
a statistical analysis of the two metrics will provide insight into the design of future EH WSNs. 

In the following figures (Figs. 6-9), we only present the numerical results for the average 
values of the two delay metrics, which have been presented in Corollaries 1-8. 

Average update age and average update cycle with different energy arrival processes: 
Figs. and 1^ show the average update age, Tua, and the average update cycle, Tuc, for different 
W, i.e., different time windows for retransmissions, W — 1, and energy arrival processes. The 
results in Figs. and are generated using Corollaries 1-4 and Corollaries 5-8, respectively. We 
can see that the different energy arrival models result in almost the same values of the average 
update age and especially the average update cycle. As the time window for retransmissions 
increases, the average update age increases monotonically and approaches its analytical upper 
bound given by Corollary 4, while the average update cycle decreases monotonically and ap¬ 
proaches its analytical lower bound given by Corollary 8. Thus, with a smaller time window for 
retransmissions, the updated status is more fresh, but the update frequency is lower. 

Average update age and average update cycle with different average harvested power: 
Fig. 8 shows the average update age, T\ja, and the average update cycle, T\jc, for different 
average harvested power values, p, with an exponential energy arrival process. The results are 
plotted using Corollaries 3 and 7. For the update age, we see that when the average harvested 
power is very low, i.e., less than —2 dBm, the update age is one time block. This is expected since 
sufficiently low average harvested power cannot enable any retransmission during time window 
W — 1, i.e., a packet is either successfully transmitted in the first transmission block right after 
the sensing block (an update age of one) or dropped due to no chance of retransmission. With an 
increase of average harvested power, retransmissions are enabled, which makes the update age 
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Fig. 6: Average update age, Tua, versus W, with differ- Fig. 7: Average update cycle, Tuc, versus W, with 
ent energy arrival processes. different energy arrival processes. 

increases beyond one. However, when the average harvested power is higher than 8 dBm, the 
average update age monotonieally deereases with an increase of the average harvested power. This 
is as expected: the sensor requires fewer energy harvesting blocks to perform retransmissions, 
and hence, the sink is likely to reeeive the paeket in a more timely manner (i.e., with a smaller 
update age). For the update eyele, we see that the average update eyele monotonieally deereases 
with average harvested power. Again, this is expected sinee a higher average harvested power 
enables more transmission bloeks within a eertain time duration, and henee, more sueeessful 
block transmissions are likely to oceur within a given time duration, i.e., the update eyele 
deereases. Also we see that when the average harvested power is very high, i.e., p > 30 dBm, 
both update age and update eyele eonverge to constant values whieh ean be obtained by letting 
p —)■ cxD in Corollaries 3 and 7, respeetively. Thus, without changing the parameters of the 
communication protocol, the improvement in delay performance is limited when increasing the 
average harvested power. 

Effect of energy cost of sensing on average update cycle: We illustrate the effeet of energy 
eost of sensing on average update cycle with exponential energy arrival process as a speeial case 
of the random energy arrival process. Fig. shows the average update age, Tua, and the average 
update eyele, Tuc, as a funetion of W, with different energy eost of sensing, Vsb- The figure 
shows that the average update age increases as W inereases (consistent with Fig. but it does 
not ehange with the energy eost of sensing, i.e., the energy cost of sensing has no effect on the 
update age. This is in perfeet agreement with our earlier observations and explanations provided 
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p, dBni 


Fig. 8; Tuc and Tua versus p with exponential energy 
arrival process. 
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Fig. 9: Tuc Tua versus W for different sensing 


power, T^sb- 


in Remark We can see that for a fixed value of W, the average update cycle increases as the 
sensing power consumption increases from 50 mW to 100 mW, i.e., the higher the energy cost 
of sensing the lower update frequency. This is in perfect agreement with our earlier observations 
and explanations provided in Remark To place these results in context with existing studies in 
the literature that commonly ignore the energy cost of sensing, we also include the result with 
zero energy cost of sensing. When Psb = 0 mW, we can see that Tuc is almost constant around 
the value of 50 and does not vary much with W. 

Tradeoff between average update age and average update cycle: Fig. shows the tradeoff 
between average update age, Tua, and average update cycle, Tuc with exponential energy arrival 
process. The different points on the same curve are achieved with different W. We can see that 
when the energy cost of sensing is comparable to or larger than the energy cost of transmission, 
e.g., Tsb = 50 mW and Tsb = 100 mW, the reduction in Tua can result in a significant increase 
in Tuc, and vice versa. For example, when Tsb = 100 mW, decreasing Tua from 15 to 5 time 
blocks, causes the Tjc to increase from 75 to 95 time blocks. However, when the energy cost 
for sensing is negligible, e.g., Tsb = 0 mW, such a tradeoff is almost barely noticeable. For 
example, decreasing Tua from 15 to 5 time blocks, results in Tuc increasing by two time blocks, 
i.e., a significant change in Tua does not result in a noticeable change in Tuc- These trends in 
Fig. 10 are in accordance with our earlier observations in Remark]^ Thus, with the consideration 
of sensing energy cost, an increase of update frequency is achieved at the expense of update 
freshness, and vice versa. 

Effect of transmit power consumption on average update age and averge update cycle. 
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Fig. 10: Tradeoff between Tuc Tua- 


Fig. lliTuo^ndTuAversusVF with different T^tx^nd^TB- 


Fig. 11 shows the impact of power consumption on Tua and Tuc, for different values of transmit 
power Ptb and RF transmit power Ptx, with an exponential energy arrival process. In reality, 
Ptb and Ptx do not have a linear relationship. Three pairs of typical values found in [25] are 
chosen. We see that both Tua and Tuc decrease with Ptb or Ptx. This is as expected: if the 
transmit power is small, Pout is high, resulting in a large number of retransmissions until the 
sensed information is successfully transmitted or IF — 1 time blocks are reached. As a result, 
Tua and Tuc are large when the transmit power is small. Thus, under these above parameter 
choices, a higher transmit power results in better delay performance. 


VII. Conclusions 

This paper has analysed the delay performance of an EH sensor network, focusing on the 
operation of a single EH sensor and its information transmission to a sink. The energy costs of 
both sensing and transmission were taken into account. Two metrics were proposed, namely the 
update age and update cycle. In order to limit the delay due to retransmissions, a time window 
for retransmissions was imposed. Using both a deterministic and a general random energy arrival 
model, the exact probability mass functions and the mean values of both metrics were derived. 
The results showed that the average update age increases while the average update cycle decreases 
with increasing retransmission window length. The average update age is independent of the 
energy cost of sensing but the average update cycle increases as the energy cost of sensing 
increases. In addition, a tradeoff between update age and update cycle was illustrated when 
the energy cost of sensing is comparable to the energy cost of transmission. Euture work can 
consider the impact of non-deterministic time for receiving the feedback signal at the sensor. 
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Appendix A: Proof of Lemma [ATI 

We first define the bloek-wise harvest-then-use proeess, and then propose and prove Lemma [ATj 
Definition A1 (Block-wise harvest-then-use process). A harvest-then-use process consists of 
energy harvesting blocks (EHBs) and energy consumption blocks (ECBs). It starts and keeps 
on harvesting energy with EHBs. Once the available, i.e., accumulated, energy is no less than a 
threshold of Q Joules, an ECB occurs, and consumes Q Joules of energy. If this condition for 
ECB is not satisfied, the process goes back to harvest energy with EHBs. 

During the harvest-then-use process, the harvested energy in the fth EHB is represented by 
^i, f = 1, 2, 3,..., and the available energy after the jth ECB is represented by Ej, j = 1, 2, 3,.... 
Due to the randomness of the energy arrival process, i.e., is a random variable, the available 
energy after each ECB, Ej, is also a random variable which only depends on ^i. Eurthermore, 
using the statistics of ^i, and modeling Ej, j = 1,2,3,..., as a random process, an important 
feature of the random process is revealed in Eemma Al. 

Lemma Al. Eor block-wise harvest-then-use process with energy threshold Q, where the 
harvested energy in each EHB, i = 1,2,3,..., are independent and identically distributed, 
each with pdf containing at least one positive right-continuous point, f{x), the available energy 
after each ECB, Ej, j = 1,2,3,...,, consists of a positive recurrent Harris chain, with unique 
steady-state distribution which is given by 

g{x) = -{l-F{x)), (A.l) 

P 

where F{x) and p are respectively, the cdf and the mean of 


Proof: The proof consists of two steps. In the first step, we prove that the energy state after 
the jth ECB, Ej, constitutes a positive recurrent Harris chain (a collection of Markov chains 
with uncountable state space). Thus, a unique steady-state distribution of Ej exists [30|. In the 
second step, we prove that ( |A.1[ ) is the unique steady-state distribution. 

Step 1: It is easy to see that the current state, Ej takes its value from a continuous state space 
and only relies on the previous energy state Sj_i, thus Ej, j = 1, 2, 3,..., forms a continuous-state 
Markov chain. Without loss of generality, we assume that sup = B, thus sup \eA<b 
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holds in this harvest-then-use proeess|^It is easy to see that the state spaee of Markov ehain Ej, 
S, is a subset of [0,5), and beeause of the harvest-then-use protoeol, any eurrent state whieh 
is higher than Q, will aeeess the interval [0, Q) in the following steps. Thus, we only need to 
prove that any state s G [ 0 , min{5, Q}) ean hit any arbitrary small interval r = in S 

with non-zero probability within finite steps. Aetually, in the following, we eomplete the proof 
with the assumption that S = [0, 5), whieh also proves that the state spaee of Markov ehain Sj 
is exaetly [0,5). 

In the following, using a eonstruetive method, we show that for Markov ehain Ej, given any 
eurrent state s G [0, min{5, Q}), there is at least a probability, q x p, that any arbitrary small 
interval r will be aeeessed with j steps, where p, q, j are defined below whieh only depends 
on the state s, the interval length r and the pdf of the harvested energy in eaeh EHB. 

Sinee pdf funetion f{x) has positive right-eontinuous points on [0,5), there exists at least 
one interval [D~,D) that satisfies [D~,D) C [0,5), D — D~ = r/2, and f{x) is positive 
right-eontinuous on [5“,5). We assume that D~ > , and the 5“ < ease ean be easily 

generated from the the following diseussions, thus is omitted due to spaee limitation. Now we 
define p = f{x)dx as the probability that harvested energy in one EHB lies in the interval 
[D~,D). Also we define f{x) = f{x) when x G [D~,D), otherwise f{x) = 0, and fi{x) is the 
i-fold eonvolution of funetion f{x). 

Thus, it is easy to see that fi{x) is positive and eontinuous in the interval {iD~,iD), and 

fi{x)dx is the probability that the harvested energy by i EHBs lies in the interval [a, 6), 
while the energy harvested by eaeh of the i EHBs lies in the interval [D~, D). Thus, letting i = 
\d{Q -f r)/r — 1] and j = [((i -f 1)5 — t+)/QJ , given the eurrent energy state s, after i EHBs, 
the accumulated energy level lies in the interval A = (s-t-z5“, s-\-iD) with positive probability 
distribution. Also we see that interval A = {jQ + t~ — D~ ,jQ + — D) C A, thus, there is 

at least (beeause we have only eonsidered the seenario that harvested energy by eaeh EHB lies 
in [5“, 5)) a probability q = inf{/- /-(a;), interval r C 5, length of r = length of A = t/2} 
that the aeeumulated energy level lies in A. Therefore, after the next EHB with probability 
p that the harvested energy lies in [5“,5), the aeeumulated energy level lies in the interval 


^Note that although we assume B is finite, the infinite case can be easily generated from the discussions below, thus is omitted 
due to space limitation. 
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[jQ + T~,jQ + r+), which means that after the eurrent state Ej = s, with j steps (eaeh step 
eonsumes the amount of energy, Q), there is at least a probability, q x p to make the Markov 
ehain hit the interval (r~, r+). Thus, Markov ehain Ej is a positive reeurrent Harris ehain [30|. 


Step 2: In the aforementioned Markov ehain, we still assume that the eurrent state Ej = s. 
Thus, in the previous state, the available energy eould be higher than Q, i.e., Ej_i = s + Q, 
and Ej_i eould also be smaller than Q, i.e., based on energy level Ej_i, there are i EHBs 
(i = 1, 2, 3,...) to make the energy level reaeh Q + s, whieh makes Ej = s. Based on the above 
and the Markovian property, the steady-state distribution of the proeess, g{x), should satisfy the 

oo _ 

following eonditions: (1) g{x) = 1 and (2) g{x) = g{x+Q) + ^ gi_i{Q+x-y)f{y)(ly. 

i=l 

where gi{x) represent the pdf of energy level after i EHBs following a ECB, whieh is given by 

g{x), i = 0 , 

(A.2) 


9i{x) = < 


gxf X f X ...X f (x), i > 0. 


V 

i convolutions 


Beeause f{x) and gi{x) > 0 for all x and 2 = 0,1,2,..., by using Tonelli’s theorem for sums 


and integrals [311, we exehange the summation and integral operator in Condition 2, thus we 
have 


Q+x 


g{x) = g{x -f Q) + 



gi{Q + x-y) f{y)dy. 


Taking ( jA.lj ) into ( |A.2| ), we have 


(A.3) 


gi{x) = - \ I F-kf ^ /... ^ / I (x) - I Fxf ^ /... ^ / I (x) I , 2 > 0. 

i—1 convolutions / \ 2 convolutions 


Sinee 0 < F(x) < 1, /(x) > 0 and /(x) = 1, when i —>■ oo, we have [30| 


(A.4) 


Fxfxf...xf (x)^0. 


i convolutions 


Erom (A.4) and (A.51, we have 


(A.5) 


1 1 

V gi(Q + X - 2/) = - (1 - F(Q + X - 2/)) + - (F(Q + x - y) - (F x f ) (Q + x - y)) 

Up p 

+ -{{Fxf){Q+x-y)-{Fxfxf){Q+x-y)) + ...= - [l-lim (f*/*/...*/) (x) | =-. 

p p \ 2^00 I ^- yj j p 

\ \ i convolutions / / 

(A. 6 ) 
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Taking ( |A.6[ ) and ( |A.1| ) into the right side of ( |A.3[ ), we have 

Q+x / ^ ^ Q+x 

^-(l-F(x + Q))+ [ -f(y)dv 

r l '' (A.7) 


g{x + Q) + 


Q+x , ^ ^ 

/ + a; - I/) j f{y)dy 


= -{l-F{x + Q)) + - {F{Q + x)- {F{x)) = - (1 - F{x)). 
P P P 


Thus, g{x) in ( |A.2| ) satisfies Condition 2. Because of — F{x)) = IE{^i} [30|, Condition 
1 is also satisfied, yielding the desired result. ■ 


Appendix B: Proof of Lemma 1 


For general random energy arrival processes, the proof is based on Lemma A1 given in 
Appendix A. First, we find an arbitrarily small Q which is a constant such that £’sb and £^tb 
are integer multiples of it. Then, from an energy perspective, we equivalently treat the proposed 
communication protocol with energy harvesting, sensing and transmission as a simple harvest- 
then-use process with EHBs and ECBs (each consumes energy, Q) as discussed in Eemma |A1[ 
Thus, the energy level after a TB, can be treated equivalently as that after a corresponding ECB. 
Therefore, the steady-state distribution of energy level after each TB is the same as that after 


each ECB, which is given in Eemma AT completing the proof. 


Appendix C: Event and Random Variable Definitions 
To assist the proofs of the main results, we use UC to denote the sequence of time blocks 


from an arbitrary STB to the next STB. Also we define two events (according to |30|) and 

several discrete random variables (r.v.s) for convenience: 

1) Event Aguc^ Given a SB, its generated information is successfully transmitted to the sink, i.e., 
STB occurs during the W blocks after the SB. 

2) Event AfaiP Given a SB, its generated information is not successfully transmitted to the sink, 
i.e., STB does not occur during the W blocks after the SB. 

3) r.v. N,1 < N <W\ Given a SB, it is followed by N TBs before the next SB. Le., if Ague 
occurs, the N TBs includes — 1 ETBs and one STB. While if Afau occurs, all the N TBs 
are ETBs. 

4) r.v. L, 1 < L <W\ After a SB, the Lth block is the last TB before the next SB. I.e., if Ague 
occurs, the Lth block is a STB, thus L is the update age. While if Afau occurs, the Lth block 
is the last ETB during the time window for retransmissions, W. 
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5) r.v. V, V > —1: Given a SB, if Ague occurs, V = —1, while if an Afau occurs, V is the 
number of the required EHBs after the time window for retransmissions, W, in order to 
harvest the amount of energy, Stb- Note that, after a Afaii, the amount of energy £^sb + 

is required to be reached in order to support the following SB and TB. Without loss of 
generality, here we assume that the energy harvesting process first meets the energy level 
Stb, and the TB consumes the energy, (V)irtually. From Lemma 1 and its proof, the 
steady-state distribution of the available energy level after the V EHBs is g(€). 

6) r.v. V, V >0. Given a SB and conditioned on a Afaii occurs, V is the number of the required 
EHBs after the time window for retransmissions, W, in order to harvest the amount of energy, 
£^tb- From the definition of V and V, it is easy to see that 


Pr{E = r;} = Pr|E = r;|Afaii} 


Pr|E = uj 
Pr {Afaii} 


,r; = 0,1,2,.... 


(C.l) 


7) r.v. E{£), E{£) > 0: Given the distribution of initial energy level, g{e), and the amount of 
target energy, 8, the required number of energy harvesting block is E{£). 

For a deterministic energy arrival process, straightforwardly we have 


Pi {E{8) =i} = l,i = 8/p. (C.2) 

For a general random energy arrival process, from the definition of E{8), Lemma and its 
proof, we have 

Pi{E{8)=%] = G,_^{8)-G,{8), z = 0,l,2,... (C.3) 


where 


G,{x) = 


1, 2 = -l, 

nx 

/ gi{u)du, i > 0, 


(C.4) 


and gi{x) is defined in ( |A.2[ ). 

For exponential energy arrival process, we know that the energy accumulation process during 


EHBs after a TB is a Poisson process [30|, thus, we have 

(8 / op e~^ ^ ^ 

Pr {E{8) =i} = - G,{8) = Pois (p 8/p) = -, 2 = 0,1,2,... (C.5) 

% • 


8) r.v. M, M > 0: Given a UC, Afau occur M times and followed by one Ague in it. 
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From the definitions of event, we know that Ague and Afaii are mutually exelusive events. Thus, 
we have 

Fsuc = Pr {Ague} and Pr {Afau} = 1 - Psuc, (C.6) 

where Ague and Afaii depends on transmit outage probability in eaeh TB, and the available 
energy after the first TB following the SB. Beeause we have assumed that the sueeess of eaeh 
transmission are independent of one another, and from Lemma 1, the distribution of the available 
energy after eaeh TB is the same, eaeh event As^c/Afaii is independent with eaeh other during 
the eommunieation proeess. Therefore, r.v. M follows the geometrie distribution 

Pr {M = m} = Psuc (1 - PsncT > m = 0,1, 2,.... (C.7) 


Appendix D: Pmf of Update Age 


From the definitions in Appendix C, the pdf of Tua can be calculated as 

= = fc = l,2,...,iy. (D.l) 

Using the law of total probability and the r.v.s defined in Appendix C, we have 

k k 

Pr{L = fc, Asuc} = '^PT{L = k,N = n,Asnc} = '^PT{N = n,E{{n - l)STB) = k- n, Ague} 

n=l n=l 

k 

= ^Pr {/7 = n, Asuc|P((n - l)£^TB) = fc - n}Pr{P((n - l)STB) = k- n} 

71=1 

k 

= 5^(1 - Pout) (Pout)""' Pr {P((n - 1)^tb) = k-n}. 


n=l 


(D.2) 


Again using the law of total probability and using ( |D.2[ ), ( |C.6[ ) becomes 

w w 

Psuc = Pr {Asuc} = ^ Pr {P = /, Asuc} = Pr {P = 1, Ague} + ^ Pr {P = /, Ague} 


1=1 
W I 


1=2 


1 - Pont + ^ ^ Pr {P = /, 27 = n, Ague} 

1=2 71=2 
VK / 

1 - P’out + X] X] Pr - 1)^^tb) = l-n,N = n, Ague} 


1=2 n=2 
W I 


1 - Pout + ^ ^ Pr {27 = n, Asuc|P((n - 1)Ptb) = / - n} Pr {P((n - 1)Ptb) = l -n} 

(D.3) 


1=2 n=2 
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W I 


1 - Pout + 5^(1 - Pont) (Pout)"-' Pr {P((n - 1)£tb) = l-n} 


1=2 n=2 


By taking ( |C.2| ), ( |C.3[ ) and ( |C.5[ ) into ( |D.2| ) and ( |D.3[ ), and then substituting ( |D.2| ) and ( |D.3| ) 
into ( |D.1| ), the pmfs of Tua for deterministic, general random and exponential energy arrival 
process are given in Theorems 1, 2 and 3, respectively. 


Appendix E: Pmf of Update Cycfe 

First, assuming that Afau occurs m times during a UC, we define r.v.s Eq, Vi, Ei,i = 1, 2,m, 
and L. Eq is the number of EHBs at the beginning of the UC until the first SB occurs which 
follows the same pmf with r.v. E{£sb + ^^tb)- Vi is the number of EHBs required to harvest the 
amount of energy, £^tb> outside the time window for retransmissions of the ith Afaii. Pj is the 
number of EHBs required to harvest the amount of energy, £^sb, following Vi EHBs after the ith 
Afaii- L is the number of blocks after a SB to the last TB before the next SB, and the TB is a 
STB. From the r.v. definitions in Appendix C, Eq, Vi and Pj follow the same distribution with 
r.v.s P(<£^sb P^^tb), V and P(£^sb), respectively, and 

Pr {l = i}= Pr {L = I, A^uc} , / = 1, 2,..., PU. (E.l) 

From Eemma 1, Eq, Vi, Ei, i = 1,2, ...,m, and L are mutually independent. 

Then, the pmf of update cycle can be calculated as 

Pr {Tuc = fc} = ^ Pr {Tuc = k,M = m} 

m 

= 'y ^ Pr •|^Po+Pi + ... + Pm + iA + --- + Pm+rn x (1 + PU) + Z/ + 1 = k,V, V, ■ ■ ■ , Vn — 

m 

rh 

= Pr |Po + Pi + ... + Pm + Ci + --- + i^ + A = ^ — rn X (1 + lU) — 1, Uf, ^mPo|, 

m=0 

k = 2,3 ,..., 
(E.2) 

where m = • For simplicity, we define the following discrete functions: 

a£,t)^FT{E{£)=i}, * = 0 , 1 , 2 , ... 
i{l) A Pr |Z = /} = Pr {L = /, Aa,J , I = 1,2, ...,W 
'(9(p) = Pr|u = u|, u = 0,l,.... 


(E.3) 
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where C(^)*) l{1) are obtained direetly from ( |C.2| ), ( |C.3[ ), ( |C.5| ) and ( |D.2[ ), respeetively, and 


'i}{v) will be derived later. Therefore, pmf of Tuc in ( |E.2[ ) ean be ealeulated as 

Pr{Tuc = ^} = X] [ C(<^SB + ^TB) fC(<^SB) * • •• - *^, *M(fc-n2(l + fE)-l). 

\ m convolutions m convolutions / 

(E.4) 

Now we derive the expression for Erom the definitions of r.v. in Appendix C, we have 

w 

i^{v) = Pr |l/ = r;} = Pr jAfau, 1/ = t;} = ^ Pr jAfau, V = v,N = n 


n=l 


= Pr Afaii, V = v,N = 1 


= Pr Afaii, V = v,N = 1 


w 

Pr I Afaii,!/ 


= V, N = n 


n=2 
M/ I 




= V, N = n, L = l 


(E.5) 


1=2 n=2 

= Pr {Afaii, iV = 1, ^(^tb) = fP + n - 1} 
w i 

+ Pr (Afaii, N = n, E{{n - 1)^tb) = I - n, E{Stj,) = W + v - 1} 

1=2 n=2 

= Pr (Afaii, N = 1\E{Stb) = hP + - 1} Pr {E{Stb) = W + v-l] 
w i 

+ Pr (Afaii, N = n\E{{n - 1)^tb) = I - n, E{Stb) = W + v - 1} x 

1=2 n=2 

Pr {E{{n - 1)^tb) = I - n, E(^tb) = fP + 

= FoatPr{^(^TB) = lE + r;-l} 

w i 

+ EE (Pout)" Pr {E{{n - 1)^tb) = / - n} Pr {P(£tb) = W + v - 1} . 

1=2 n=2 

By taking funetions ( |E.5[ ), and i{l) in ( |E.3[ ), into ( |E.4[ ), and letting ( |C.2| ) and ( |C.3| ) 

substitute Pr{P(£^) = z}, the pmf of Tuc for deterministie and general random energy arrival 
proeess ean be ealeulated, respeetively, as given in Theorems 4 and 5. While for the exponential 
energy arrival proeess, by using the sum property of Poisson distribution, we have 

Pr {P(^i)f + E{£ 2)2 = *} = Pr {P(^i + £ 2 ) = 1 } , (E.6) 

where E{£i)i and ^(£’ 2)2 are two independent random variables whieh have the same distribution 


with E{£i) and ^(£^ 2 ) defined in Appendix C, respeetively. Therefore, letting (|C.5|) substitute 
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Vi{E{£) = i}, the pmf of Tuc for exponential energy arrival proeess ean be further simplified 
as given in Theorem 6. 


Appendix F: Average Update Cycle 
Based on Appendix E, average update eyele ean be ealeulated as 

OO 

Tuc = E {E {Tuc|M}} = ^ Pr {M = m} E {Tuc|M = m} 

m=0 

OO 

= ^ ^ Pr {iff = 77i} E {Uo+-Ei + - ■ ■+Um + Ui + V 2 + ' ■ ■ + 14ra+??r x (1 + 1U) + 1+Tua} 

m=0 

OO 

= ^ ^ Pr {iff = 77i} (E {Uo}+E {Ui} + - ■ 'AE {Uml+m x V -\-m x (W + 1)+7ua + 1) • 

m=0 

(FT) 


From Appendix E, we have 


TIT r 1 ^SB + £tb in, r to 1 ^SB . . „ 

E{Uo| =- - -, E{Ui| = —,^ = 


(F.2) 


After taking ( |E.5[ ) and ( |C.6| ) into ( |C.1[ ) and some simplifieations, the expeetation of V ean be 
ealeulated as 


1 _ p 

-L c;i 




out 




W-2 


W-2 


TB 


v =0 


i-P«, 


p 


5^zPr{U(^TB) =0-(W^-l) l-5^Pr{i?(^TB)=*} 


i=0 


i=0 


W I 


1-Ps, 


5^ (Pout)” Pr {E{{n - 1)Stb) = I - n} x 


1=2 n =2 


& 


W-l-1 


W-l-1 


TB 


P 


tPT{E{STB) = t}-{W-l)[l- Pr{P(^TB)=*} 


i=0 


i=0 


(F.3) 


By taking ( |E2| ), ( |E.3[ ) and ( |C.7[ ) into ( |E.1[ ), and further substituting Pguc and Tua given in 
Corollaries 1, 2 and 3, average update eyele for deterministie, general random and exponential 
energy arrival proeesses are given in Corollaries 5, 6 and 7, respeetively. 


Appendix G: Asymptotic Eower/Upper Bounds 

Erom Corollaries 1 and 2, it is easy to see that Tua inerease with W. While for Tuc, the 
monotonieity is not explieitly observed from Corollary 6. Due to spaee limitations, a sketeh of 
the proof is given: When W inereases, more TBs are allowed, thus more STBs oeeurs during 
the eommunieation proeess, whieh also means shorter average update eyele. 
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When W ^ oo, the sensed information in each SB will be successfully transmitted to the 
sink, i.e., Ague always occurs and Pguc —^ 1- Thus, UC contains the EHBs to harvest the amount 
of energy, £^sb + the SB, and the blocks in Tua- Based on this explanation, for the average 
update age, we have 


lim fuA = VPr{iV = n}E{TuA|iV = n} = V(l-Pout)(Poutr“'E{n + P((n-l)^TB)} 

W^OO ^ ^ 

n=l n=l 

= f;(l - P„„,) (A,,.)- (n + (n - 1)^) = 1 + + li 

n=l ^ P / \ P / 


1-P 


For the average update cycle, we have 

oo 

lim Tuc = Pr {N = nj E {Tuc|A^ = n} 

W^oo 


(G.l) 


n=l 


- Pout) (Pout)"-' E {P(^SB + ^tb) + 1 + n + P((n - 1 )^tb)} 


n=l 

oo 


y~^(l~-Pout) (-Pout) ( ^^+14 


n=l 


^^SB + uStb 
p 


= 2 + 


^TB Pout ^SB + ^TB 


P 


1-P 


out 


p 


(G.2) 
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